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When two fully polarized ferromagnetic (F) wires with opposite polarizations make con- 
tact with a spin-singlet superconductor, a potential-induced current in wire 1 induces 
a non-local current of equal magnitude and sign in wire 2. The magnitude of this cur- 
rent has been studied in the tunneling limit and found to decay exponentially with the 
distance between the contact. In this paper we propose a new structure in which this 
novel non-local effect is increased by orders of magnitude. We study the spin-dependent 
electronic transport of a diffusive nanoj unction and demonstrate that when a normal 
diffusive region is placed between the F leads and superconductor, the non-local initially 
increases with the separation between the F leads, achieving a maximum and decays as 
a power law with increasing separation. 



Experimental studies of electronic transport properties of nanostructures con- 
taining both ferromagnets (F) and superconductors (S) 1 . 2 . 3 ^ 4 ! 5 ^,7,8,9,io,n_ revea [ 
novel features, not present in normal-metal/superconductor (N/S) junctions, due 
to the suppression of electron-hole correlations in the ferromagnet. When spin- flip 
processes are absent, further effects are predicted, including the suppression of the 
conventional giant magnetoresistance ratio in diffusive magnetic multilayers and 
the appearance of non-local currents when two fully polarized ferromagnetic (F) 
wires with opposite polarizations make contact with a spin-singlet superconductor 
13,14 . The latter effect has been highlighted, because of interest in the possibility of 
generating entangled pairs of electrons at a superconductor (S) interface 15 ' 16 - 17 . A 
recent study of such a junction in the tunneling limit predicts that the magnitude 
of the non-local current will decrease exponentially with the distance between the 
F contacts. It is therefore of interest to ask how the magnitude of this novel effect 
can be enhanced. 

In this paper we propose a hybrid nanostructure in which the non-local current 
is enhanced by orders of magnitude compared with the geometry of ref. 13 . The 
proposed structure is sketched in figure 1 and comprises two clean F wires, each of 
width Mf, separated by a distance M, in contact with a diffusive, normal metallic 
region of area A = L.(2Mf + M), which in turn makes contact with a spin-singlet 
superconductor . 

In the linear response limit, the current- voltage relation of a such a hybrid 
structure connected by non-superconducting wires to normal reservoirs was first 
presented in ref 18 . If ij.,/2 are the currents leaving reservoirs 1 and 2, and vi,V2 
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their respective voltages, then at zero temperature 18 

h = [2e 2 /h][{N -R + R a )( Vl - v) + (T' a - T>)(v 2 - v)] (la) 

and 

h = [2e 2 /h][(T a - T )(vi — v) + (N — R + R' a )(v 2 - «)], (16) 

where v is the condensate potential. In this expression Ro is the coefficient for 
an electron from reservoir 1 to be reflected as an electron back into reservoir 1, 
while Tq is the coefficient for an electron from reservoir 1 to be transmitted as an 
electron into reservoir 2. R a is the coefficient for an electron from reservoir 1 to be 
Andreev reflected as a hole into wire 1 and T a is the coefficient for an electron from 
reservoir 1 to be transmitted as a hole into wire 2. Finally R' , Tq, R' a , T' a are the 
corresponding coefficients for electrons originating from reservoir 2. All coefficients 
are evaluated at the Fermi energy and are obtained by summing over the elements 
of sub-matrices of the quantum mechanical scattering matrix S. As a consequence 
of unitarity of S, they satisfy N = R Q + R a + T a + T = R' + R' a + T' a + T{>, where 
N is the number of open scattering channels in the wires. 

As discussed in ref 19 , equations (1) are a convenient starting point for describing 
a wide variety of transport phenomena in hybrid superconducting nanostructures. 
For example if lead 2 is employed as a voltage probe, with I 2 set to zero, then 
equation (lb) yields 

(v 2 - v)/( Vl — v) = (N — R' + R' a )/(T - T a ). (2) 

In the absence of Andreev processes, this ratio is positive. However if Andreev 
transmission dominates normal transmission, the ratio can be negative, leading to 
novel negative 4-probe conductances in hybrid superconducting structures 20 . The 
appearance of non-local currents is similarly obtained from equations (1) by setting 
v 2 — v = 0, which yields 

h = [2e 2 /h] [N-R + R a ] ( Vl - v) = [2e 2 /h] [2R a + T + T a ] («i - v) (3) 

and 

I 2 = [2e 2 /h][T a -T ](v 1 -v). (4) 

When all Andreev processes are absent, this yields the expected result I 2 = —h = 
— (2e 2 /h)To(vi — v), whereas when Andreev reflection and normal transmission are 
completely suppressed, I 2 = +h = (2e 2 /h)T a (v\ — v). The latter occurs when the F 
wires are completely spin polarized and the polarization in wire 1 is the opposite of 
that in wire 2. More generally I\ and 1 2 will have opposite signs whenever Andreev 
transmission dominates normal transmission. As noted in ref , this can occur even 
if the current carrying wires are not ferromagnetic. 

In the absence of a normal region in front of the superconductor (ie L = 0) 
and in the presence of tunnel junctions at the F-S interfaces, the analysis of ref 13 
predicts that both I\ and I 2 decay as exp —(2L/tt£,), where £ is the superconducting 
coherence length. To analyze this structure for finite L we solve the Bogoliubov - 
de Gennes equation on a square tight binding lattice and compute the scattering 
matrix using an exact recursive Green's function technique. 
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We start by defining a tight-binding lattice of sites with the geometry of figure 
1. Each site is labeled by a lattice vector / and possesses particle (hole) degrees 
of freedom ip a (l) (^(O); with spin a = ±1. Ferromagnetism is incorporated via a 
Stoner model, with an exchange splitting ho(l) on site I. In the presence of local s- 
wave pairing described by a superconducting order parameter A(Z), the Bogoliubov 
equation takes the form 

(0 - \4) - oh (l) - J2 TV (*"*+ + ^A(f)0- ff (F) (5) 

6 

E^(l) = -[e(l) - ah (l)W(l) + E^( f + % - ^A*(f)V>- CT (f), (6) 

where I + S sums over the nearest neighbors of I. 

If I belongs to the normal diffusive region of figure 1, e(l) is chosen to be a 
random number, uniformly distributed over the interval e — W/2 to e + W/2, 
whereas in the clean F and S regions = e . In the S region, the order parameter 
is set to a constant, A(Z) = Ao, while in all other regions, A(F) = 0. If Fis located 
in one of the clean ferromagnetic leads, then Stoner splitting is h (l) — hj, where 
j = 1,2 labels the lead containing site I. For aligned moments h\ = h 2 = cm and for 
anti-aligned moments, hi = —h 2 = cm- If I is not located in one of the ferromagnetic 
leads then ho(l) = 0. The nearest neighbor hopping element 7 fixes the energy scale 
(ie the band-width), whereas eo determines the band-filling. In what follows, to 
model the experimentally-relevant regime of Ao <C Fermi energy(= 4-y — eo) and 
Mf > Fermi wavelength, we set 7 = 1, e = 0.2, e M = 3.8, M f = 10 and A = 0.1. 

By numerically solving the Bogoliubov - de Gennes equation for a given energy 
E, the scattering matrix S can be computed for a given realization of the disorder 
and average values of measurable quantities obtained ensemble averaging the quan- 
tity of interest. For a given value of disorder W, the elastic mean free path A e ; of the 
normal region is obtained by from a separate calculation of the ensemble-averaged 
dc conductivity of a rectangular normal bar, connected to normal reservoirs. A e ; is 
then obtained by comparing the result with the Drude formula. In what follows, 
the lattice constant is set to unity and W chosen such that A e ; = 10. 

For E — 0, figures 2 and 3 show the dependence on L and M of the ensemble- 
averaged Andreev transmission coefficient in the presence of non-magnetic wires, 
while figures 4 and 5 show the corresponding quantities in the presence of fully 
spin-polarized wires with opposite magnetizations. For small L (eg L = 1) and 
finite M, figures 3 and 5 show that T a is exponentially small, as predicted by 13 . 
However as L is increased to a finite value, T a shows a dramatic enhancement. 

Although this structure is difficult to describe analytically, the qualitative be- 
haviour of T a is clear and can be understood in terms of a simple resistor model, in 
which T a is proportional to the conductance of electrons from lead 1 travelling first 
to the N-S interface and then to lead 2. Consider for example figure 2, which for 
M 7^ shows the expected exponential suppression of T a at L = 0. For small L, T a 
increases linearly with L. The slope of this linear regime is independent of M for 
small M < X e i, reflecting the ballistic nature of the N region in this limit, whereas 
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for the opposite limit of M > A e ;, the slope decreases with increasing M, reflecting 
the diffusive behaviour of the N region. In the diffusive regime, for M << L, this 
leads us to expect the T a is proportional to L/(M + 2M/), whereas for L << M 
it is proportional to (M + 2Mf)/L, which suggests that the maximum value of T a 
occurs when L w M + 2Mf. Figure 4 shows that same qualitative behaviour also 
persists in the presence of ferromagnetic leads. 

In conclusion, we have shown that in contrast with the tunneling regime, where 
T a is exponentially small, the non-local current can be dramatically enhanced by 
including a normal region between the F leads and superconductor. In this case the 
non-local current no longer decays monotonically with M or L, but instead exhibits 
a maximum corresponding to a square normal region with L rs M + 2Mf. 
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Figure 1: The hybrid nanostucture analysed. 
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Figure 2: For normal leads (ie when em = 0) this figure shows Andreev transmission against 
length L for different contact separations M . 
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Figure 3: For normal leads, this figure shows Andreev transmission against contact separation M 

for different lengths L. 




Figure 4: For ferromagnetic leads (ie cm = 3.8J with anit-aligned magnetiszations, this figure 
shows Andreev transmission length L for different contact separations M . 
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Figure 5: For ferromagnetic leads (ie €m = 3.8) with anit- aligned magnetizations, this figure 
shows Andreev transmission against contact separation M for different lengths L. 
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